Proposal for an Experiment to Test a Theory of High 
Temperature Superconductors 

CM. Varma 
Bell Laboratories, Lucent Technologies 
Q\~- Murray Hill, NJ 07974 

On 



O 

o 



o 



c3 



o 



3 



Abstract 



A theory for the phenomena observed in Copper-Oxide based high tempera- 



J> ■ ture superconducting materials derives an elusive time-reversal and rotational 

o 



symmetry breaking order parameter for the observed pseudogap phase ending 
at a quantum-critical point near the composition for the highest T c . An ex- 



0^ ■ periment is proposed to observe such a symmetry breaking. It is shown that 



Angle-resolved Photoemission yields a current density which is different for 



i-rt ■ left and right circularly polarized photons. The magnitude of the effect and its 

a: 

momentum dependence is estimated. Barring the presence of domains of the 



predicted phase an asymmetry of about 0.1 is predicted at low temperatures 
in moderately underdoped samples. 



I. INTRODUCTION 

Despite twelve years of intesive experimental and theoretical studies of copper-oxide 
based superconducting compounds, |l[ no consensus about the fundamental physics or even 
about the minimum necessary Hamiltonian to describe the phenomena has emerged. One 
of the few theoretical ideas which has clearly survived experimental tests is that at density 
x ~ x c near that for the maximum superconducting transition temperature, the normal 
state is a marginal Fermi-liquid (MFL) 0. The MFL is characterized by a scale- invariant 
particle- hole fluctuation spectrum which is only very weakly momentum dependent. One of 
the predictions of MFL hypothesis is that the single-particle spectral function G(k, uj) has 
a nearly momentum-independent self-energy proportional to max(\uj\,T). The frequency 
and temperature dependence as well as the momentum independence have been tested in 
angle-resolved photoemission experiments. 

The observed non- Fermi-liquid behavior near x ~ x c in resistivity, thermal conductivity, 
optical conductivity, Raman scattering, tunneling spectra, and the Cu nuclear relaxation 
rate follow from the MFL hypothesis. The scale-invariance of the MFL fluctuations implies 
that a quantum-critical-point (QCP) exists at x = x c , near the optimum composition. One 
expects that, in two or three dimensions, the QCP at T = is the end-point of a phase of 
reduced symmetry as x is varied. Similarly a line of transitions or at least a cross-over is 
expected at a finite temperature T p (x) terminating at (x = x c ,T = 0). Indeed the generic 
phase-diagram, Fig. (1), of the copper-oxide compounds around x ~ x c displays such a 
topology. Region I has MFL properties dominated by quantum-fluctuations, Region III 
displays properties characteristic of a Fermi-liquid, while Region III - the pseudo-gap region 
displays a loss of low-energy excitations compared to Region II. Below the line T p (x) between 
regions I and II, the single-particle spectrum displays lowered rotational symmetry, while 
no translational symmetry appears broken. The superconductivity region sits spanning the 
three-distinct normal state regions. 

Fig. (1) may be compared to the topologically similar phase-diagram of some heavy- 



Fermion compounds, in which the line T p (x) corresponds to an antiferromagnetic transition. 
H From this point of view the crucial question in Cu-0 compounds is the symmetery in 
Region II of the so-called pseudo-gap phase. 

A systematic theory |||| starting with a general model for Cu-0 compounds provides 
an answer to this question. The region (II) in Fig. (1) is derived to be a phase in which a 
four- fold pattern of current flows in the ground state in each unit cell as shown in Fig. (2). 
Time-reversal symmetry as well as rotational symmetry is broken but the product of the 
two is conserved. This phase has been called the circulating-current (CC) phase. Quantum 
fluctuations about this phase are shown to have MFL fluctuations, characteristic of Region 
I. The same fluctuations promote "d" or generalized "s"-state pairing depending on the 
Fermi-surface at a given doping. 

While a microscopic theory in agreement with most of the principal experimental results 
has been presented, one can be confident of the applicability of the theory only if the 
CC phase is directly observed. The CC phase has a very elusive-order parameter. The 
four-fold pattern of microscopic magnetic moments in each unit cell changes the Bragg 
intensity for polarized neutrons at certain pre-existing Bragg spots. But the intensity for 
nuclear-scattering at these Bragg spots is O(10 4 ) the predicted magnetic intensity. Muon 
spin-resonance (/x-SR) would be a possible probe, but the magnetic field from the current 
pattern in Fig. (2) is zero at most symmetry points and along the principal symmetry lines, 
where muons are known to sit preferentially. Perhaps, an additional perturbation such as an 
external magnetic field can be used to lower the symmetry at the sites preferred by muons. 
In that case \i — SR could be used to search for the predicted phase. 

I propose here a new kind of experiment, which is a microscopic analog of circular 
dichorism. The idea is that ARPES at a specific point near the Fermi-surface should have 
an electron yield which is different for right circularly polarized and left circularly polarized 
photons if the ground state has T-breaking of the form shown in Fig. (2). Further the 
relative intensity should change in a systematic fashion with the momentum around the 
Fermi-surface. 



I present below the results of the calculation based on this idea and then discuss the 
feasibility of the experiment. The idea itself is more general than the specific application 
to copper-oxides. Any time reversal breaking phase will in general yield a different current 
density for right and left circularly polarized photons. (But the characterestic signature 
of the state, revealed by the momentum dependence of the asymmetry in the current for 
the left and right circular polarizations, must be calculated anew for each possibility.) The 
experiment may for example be tried to see if the superconducting state of the compound 
Sr2Ru04, breaks time-reversal symmetry. 

II. ARPES WITH POLARIZED PHOTONS 

My object is to deduce the polarization and symmetry dependence of ARPES current and 
a rough estimate of its magnitude. For this purpose, a simple calculation using tight-binding 
wave-functions in the solid is sufficient. 

Assume that a photon of energy uj shone on the crystal produces a free-electron with 
momentum p and energy E p at the detector due to absorption of the photon by an electronic 
state |k > inside the crystal of the crystal momentum k. The momentum of the photon is 
assumed very small compared to k and p. The current J pk collected at the detector for 
uniform illumination over a given area is || 

J p , k = 2tt e f (e k ) | ( p \H'\ k) | 2 5 (E p - e k + hu) (1) 

where /(e k ) is the Fermi-function. 

The primary contribution of the current is from the matrix element 

(p\H'\k) = — /dre ipr A-W k (r) (2) 

2mc J 

where A is the vector potential of the incident photons and \I/ k (r) is the wave function of 
the state |k >. There is a smaller contribution due to the gradient of the potential at the 
surface which is briefly discussed at the end. 



A. Wavefunctions 

The creation operator for the tight-binding wavefunctions for the conduction band of 
Cu-0 metals (assumed to be a two-dimensional metal) for the case that the difference in 
energy of the Cu-d x 2_ y 2 level e^ and the 0-p X}V levels e p is much less than their hybridization 
energy, and when the direct Oxygen- Oxygen hopping parameter t pp is set to zero are 

' k ^S + < f 4£^) (3) 

where s XtV = smk x>y a/2 and s\ y = sin ^ xCi + sin 2 ya . Spin labels have been suppressed. 
d£, vtx, y are respectively the creation operators for the basis wave-functions 

^(k) = ^£ e " ik ' Rl ^( r - R 0, 



v^) = "4 E e- ?k - Ri e-^H P ^{v - Ri - ^) (4) 



'N 

where 0d(r — Ri) is the d x 2_ y 2 atomic orbital at the Cu-site Ri and 4> px (r — Rj — ^j is the 
p x wavefunction at the oxygen site at Rj + ^f, etc. 

In the circulating current phase, the conduction band wave-function is modified to [EJ 



|k) = (|k o ) + e |k 1 ))/ A /l + ^ s 2 s 2 (5) 

where 

|ki) ~ s x s y (s y vt x - s x pt y ) /s xy . (6) 

In Eq. (5) 9q characterises the strength of the symmetry-breaking. 

B. Matrix-elements and Current 

In order to evaluate the matrix element, Eq. (2), I write 

Mr) =^(|r|)^A 
PM ( r ) =M\ r \)fa V = x,y 



(7) 



^(i(l r l) and ^p M (|r|) are characterized by a fall-off distance a of the order of the atomic size. 
Then 



v<Ur) - M 2t fiFi& + £f ^(M) 



V0 d (r) » f| (^±M±££) (^ + ^Mll ^(| r |). 



(8) 



|r| ) \r\ \r\ 

-yy+zz \ (x 2 - 
|r| J r 

We need the (two-dimensional) momentum distribution of the wavefunctions in Eq. (2). 
For this purpose define 

dre *(p.^»)V^(r) = i fZ(p x ,py); v = x,y . (9) 

Note that fd(.Px,P y ) can be written as the product of an odd function of p x and an even 
function of p y , etc. Similarly, 

J dv e^ x+ ^V A P M = flippy) (10) 

fp {Px,Py) is the product of an even function of p x and an even function of p y , whereas 
fp (PxiPy) i s the product of an odd function of p x and an odd function of p y . The definitions 
in Eqs.(9,10) ensure that all the f's are real. The f(p)'s fall off for p of order the inverse 
atomic size. Therefore for p's in the first or second Brillouin zone they are approximately 
constant. 

In terms of these quantities, the matrix element in Eq. (2) is calculated. Consider 
the case of left and right circularly polarized photons with vector potentials A^ and A r 
respectively 

A v = A(x±iy) . (11) 

Then a straight forward calculation leads, to leading order in # to 

<P 1^1 k >^ = ( *fc) A ^Gx, Gy S(p - k - G) {(Ro(G, P , k) ± iI (G, p, k)) (12) 

+6 (±i?!(G, p, k) + i h(G, p, k))} (13) 

where G = (G x ,G y ) are the reciprocal vectors, and 



R (G, p, k) = /*(p) + (g(G y , k x )f» x (p) + g(G y , fc x )#,(p)) (14) 

/ (G,p,k) = /l(p) - (g(G x ,k y )Jl(p)+g(G y ,k x )fl(p)) (15) 

J*i(G,p,k) = (sin 2 (^ g(G x ,ky)fc(p) -sin 2 ^ (G y A)/* (p)) (16) 

/i(G,p,k,k) = (sm 2 (^fj g(G x ,k y )fy px (p) -sin 2 ^ <KG S , *»)#,) (17) 

In the above 

/ x sin(ra/2) 

g(r, s) = -, (18) 

J sin 2 (ra/2) + sin 2 (so/2) 

In the usual experimental geometry, the contribution from each G is selected separately. 
For a particular G, the current with polarization £ or r to first order in 9 is 

J/, r (G,p) = -^ fte + / c 2 ) ±29(R R 1 + VOl ( 19 ) 

so the relative asymmetry of the current, 

S(G,p) = (J < -J r )/^ + J r ), « 8^ (i2 i2 1 + / / 1 )/(^ + / 2 ) (20) 

III. DISCUSSION OF ARPES - ASYMMETRY 

Eqs. (19) and (20) are the principal result of the calculation. It is worthwhile noting 
several aspects of the predictions. For G x = G y = 0, when only the d-orbitals contribute 
to the photo-current, S = for all k. For G x = G y , the asymmetry vanishes along the 
zone-diagonal k x = k y and is maximum for the zone-boundaries (k x a = n, k y a = 0); (k x a = 
0, k y a = it) with a smooth variation in between. Asymmetry patterns for other G's may be 
obtained from Eqs. (14)- (17). 

In Ref. (6), 9q is estimated to be O(10 -1 )(a; c — x) 1 ! 2 for x < x c . So at x c — x ~ 5 x 10~ 2 , 
the asymmetry is predicted to be O(10 _1 ), at T ss T p (x). 

7 



The proposed experiment is to measure the ARPES current in under-doped samples for 
a fixed relative geometry of the incident photon-beam, crystalline surface, and the detector 
to select a p and G and then simply switch the polarization of the incident photons, and 
measure the current again. The experiment should then be repeated for different p and G. 
The effect should set in for temperatures T < T p (x) and have a momentum dependence 
predicted by Eq. (20) and Eqs. (14)-(17). 

The principal difficulty of the experiment is the possible presence of domains of the CC 
phase. The domains consist of regions in which (#o) in the wave-function (5) is replaced by 
(— #o)- This leads to a mutual switching of the pattern of the current within the unit cells 
(and a current flow along the domain boundary). The effect calculated in Eq. (20) to linear 
order in 9$ then averages to zero for equal number of the two-kinds of domains in the surface 
area S from which the current is collected. If the characteristic domain size is D, an effect 
proportional to (D/S) 1 ^ 2 6q is still to be expected. Also, Eqn. (1) yields asymmetry terms 
proportional to 6q, which are not affected by the domains. However, these may be too small 
to be observable. 

In the above, circularly polarized photons, with the plane of polarization along the 
surface of the crystal have been considered. There is also an effect linear in 6 for photons 
linearly polarized along the normal to the surface due to the potential gradient at the surface 
(VV) S . This effect, proportional to (VV) 2 changes sign for a given G x = G y , as p x and p y 
are interachanged in a d-wave like fashion. Observation of this effect requires rotating the 
sample. It is also affected by domains. 
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FIGURE CAPTIONS 



Fig. (1): Generic phase-diagram of the cuprates for hole-doping. 

Fig. (2): Current pattern predicted in phase II of Fig. (1) in References (5) and (6). 
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